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On Sensitivity of Cherenkov Radiation to the Dynamics of High Energy 
Cosmic Ray Interactions 
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A simulation study of the effects resulting from creation of quark-gluon plasma in high energy collisions of 
cosmic iron nuclei with the air has been carried out. The Cherenkov light emission has been found to reveal some 
features of the dynamics of high energy hadronic interactions. 



1. Introduction 

One possible way to investigate rare ultra- 
energetic collisions of cosmic ray primaries with 
the air is to investigate the Cherenkov light emit- 
ted in atmospheric showers. Charged daugh- 
ter particles produced along a hadronic cascade 
have mostly very high velocities and therefore 
emit in the air Cherenkov photons which in prin- 
ciple carry an information about longitudinal 
shower development including an information on 
the characteristics of the primary collision, e.g. 
particle multiplicity. As optical photons suffer 
relatively smaller absorption in the atmosphere, 
their density on the ground is higher than that 
of charged particles and their lateral distribution 
is broader. The Cherenkov light highlights dif- 
ferences in the mechanism of hadronic collisions 
that are expected to occur in the first stage of 
the shower development and that can propagate 
further during development of the cascade. 

2. Primary Interaction 

Since available experimental data do not con- 
tradict to the fact that the abundance of 56 Fe 
nuclei in the primary flux may be above 10% pQ, 
a non-negligible portion of iron-nitrogen or iron- 
oxygen nucleus-nucleus collisions can be observed 
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in upper atmosphere. The number of nucleons 
participating in such collisions may well grow up 
to a few tens and a fireball of hot and dense ele- 
mentary matter of quarks and gluons (QGP) can 
be created. Its thermal equilibrium is achieved in 
rapid statistical energy distribution among par- 
tons. As a result an almost ideally thermal- 
ized system of nucleons, pions and kaons (heav- 
ier baryons and mesons are neglected in our ap- 
proach) is produced in the primary collision. 

The energy density of the fireball can be esti- 
mated from the Bj0rken formula jHJ, e = ^rpr 
where R is the radius of the smaller nucleus, t£e 
formation time was chosen as r ~ 1 fm/c and 
AE is c.m.s. energy of produced particles in the 
unit rapidity interval in the vicinity of its c.m.s. 
value. Energy densities and corresponding upper 
bounds of the transition temperatures achieved 
in our simulations of the primary interaction of 
the cosmic iron nucleus with the air nitrogen are 
shown in Fig^ as functions of the Fermi energy 
variable F w (satat)*, where ^/snn is an invari- 
ant mass of a pair of interacting nucleons [2]. 

When the number of wounded nucleons -/Vr nt in 
the first collision is greater than a certain value 
(typically N\ nt > 20), we assume that several 
nucleons melt down into the QGP phase and 
a certain part of the total mass of the fireball 
Ehot ~ V^/sf is released. The model parameter rj 
(rj = 0.1 — 0.4) is assumed to be independent of 
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the collision energy and the system size. A por- 
tion of released energy is channeled to nucleons 
with momenta generated according to the rela- 
tivistic Boltzmann distribution isotropically with 
respect to the fireball center of mass system. The 
remaining energy is then used to create a number 
of 7r and K mesons with isotropic and relativistic 
Boltzmann distributed momenta in fireball sys- 
tem. Thermalized hadronic species are produced 
with probabilities N : it : K « 10 : 85 : 5 and 
with approximately equal abundances of different 
charges. Chemical freeze-out is controlled by a 
constant temperature (kT — 170, 340 MeV) and 
zero chemical potentials. Multiplicities of light 
mesons produced in the initial interaction, as cal- 
culated in our model, are shown in FigEI 

3. Shower Simulations 

To achieve the quantitative results, the COR- 
SIKA simulation code 4 has been adapted to 
calculate the contributions of the Cherenkov pho- 
ton intensity stemming from the quark-gluon 
equilibration in primary collisions of the cosmic 
iron nuclei with the air. The QGSJET01 model 
of high energy interactions has been employed; 
the GEISHA procedure has been used to treat 
hadronic collisions at energies E < 80 GeV. It is 
worth to note that the uncertainty related to the 
modelling of standard hadronic interactions turns 
out to be of minor influence, as the present inves- 
tigation of QGP effects relies upon comparative 
characteristics of the detectable photons. 

We analyze the contributions of Cherenkov 
light to the lateral ground level distribution func- 
tions of iron-induced showers for a given primary 
energy (Ep c — 100 TeV — 10 EeV) averaged over 
100 showers. The primary zenith angle was fixed 
at zero degrees; the altitude of the initial interac- 
tion was left free. The kinetic energy cutoffs for 
hadrons and muons have been set to 0.3 GeV; for 
electrons/positrons and photons we use cutoffs of 
20 MeV and 2 MeV, respectively. In all calcula- 
tions a thinning level 10~ 7 has been adopted. The 
Cherenkov light has been collected by a row of 100 
ideal detectors with detection area 3x3 m 2 and 
spacing 10 m located in north-south direction at 
a height of 1 km above the sea level (813 g-cm~ 2 ). 



In all simulations the standard U.S. atmosphere 
has been used; propagation of the Cherenkov light 
has been corrected for attenuation in the air. 
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Figure 1. Estimated energy densities and transi- 
tion temperatures achieved at the first Fe-N in- 
teraction are depicted as functions of the Fermi 
energy variable (see text). 
A 

V 



A 

V 
60 

3 



▲ it kT=3'10 MeV 
to kT=1 70 MeV 



77 = 0.4 



▼ 

A 



AT 



<7T>, 

<K>, 



<N> 
<N> 



Log F(GeV m ) 

Figure 2. The average yields of pions and kaons 
per average number of participating nucleons as 
functions of the Fermi energy variable. SPS and 
RHIC experimental results for pions are indicated 
by the straight line. 

4. Effects of QGP formation 

In our simple model, more secondary particles 
are created in the first interaction when QGP is 
formed (QGP shower) compared to the ordinary 
Fe-N collisions (Fe shower). Air showers follow- 
ing QGP production have their maxima higher 
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in the atmosphere and produce noticeably more 
Cherenkov light at higher altitudes; this effect 
grows with increasing incident energy. In Fig|3] 
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Figure 3. Average ground level densities and rel- 
ative QGP/Fc numbers of Cherenkov photons are 
depicted as functions of the slant depth of their 
production. Hatched (dotted) boxes correspond 
to maxima of the QGP (Fe) showers. 
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Figure 4. Relative QGP/Fe ground level lateral 
distribution of Cherenkov photons. Average val- 
ues of QGP and Fe shower maxima are indicated, 
the average ground level densities and ratios of 
total numbers of Cherenkov photons produced 
along the QGP and Fe showers are depicted as 
functions of the slant depth of their birth. 



Average ratios of the ground level Cherenkov 
photon lateral densities from the QGP and Fe air 
cascades are depicted in FigQJ The photon core 
densities of the QGP showers are typically by a 
factor of 1.2 smaller than the core densities of the 
Fe showers of the same primary energy. At very 
high primary energies, Ef c > 1 EeV, the relative 
fraction of photons at outer parts of the QGP 
shower increases. Thus the QGP effects could 
be visible comparing the number of Cherenkov 
photons reaching the detector level relatively far 
from the shower core with all detected photons. 

Any influence of the transverse momenta of the 
direct hadrons in the first interaction becomes 
negligible as in the case of very energetic primary 
the lateral profile of shower is dominated by inter- 
actions of particles at lower energies in the later 
part of the shower development. The QGP show- 
ers behave similarly as ordinary hadronic showers 
initiated by the iron primary of somewhat lower 
energy. However, the lateral spread of Cherenkov 
photons caused by the equilibrated QGP prod- 
ucts is larger. All these effects depend on both the 
primary energy and the distance between detec- 
tor and shower maximum; the latter observable 
is also energy dependent. 

5. Conclusions 

The lateral ground level Cherenkov photon 
densities of the extended air showers can be used 
to identify showers that follow the QGP transfor- 
mation in the first interaction. It is worth noting 
that similar QGP related effects can be observed 
in lateral distributions of high momenta muons 
detected underground [S]- 
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